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ABSTRACT 


Laminated composite plates have been analysed for their 
vibrational characteristics using a new technique based on the 
Gram-Schmidt procedure, to generate the orthonormal bases for 
the series expansion of the displacements, in the Rayle igh-Ri tz 
method. The technique has been demonstrated to be very powerful- 
for the analysis of moderately thick plates even shear 
deformation is neglected. The rapid convergence, the ease with 


which starting terms 

maybe chosen and 

the simple 

Gaussian 

integration possible. 

make the 

present 

technique 

extremely 

useful. Finally, the 

technique has 

been 

effectively 

used to 

analyse free vibrations of tapered 

plates . 
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CHAPTER 1 

INTRODUCTION AND LITERATURE SURVEY 


£1.13 Introduction: 

The primary load carrying and supporting framework of any 
system, be it a building, an automobile or an aircraft, is its 
structure. In particular. Aerospace structures have always 
afforded the greatest of challenges to the designer due to the 
innumerable constraints he has to work with. Any Aerospace 
structural design has to satisfy - apart from its main 
functional requirement - several other criteria such as maximum 
strength to weight and stiffness to weight ratios’, high 
fatigue strength, and good reliability. 

In order to be able to predict the total response of a 
structure to any input, a knowledge of the responses of 
individual structural elements that go to make up the overall 
structure, is essential. Some typical structural elements are 
solid and thin-walled beams, plates, shells, struts and the 
like. The response of these elements with varying degrees of 
complexity of material, loads and boundary conditions needs to 
be fully understood before an actual design can be attempted. 

The advent of Composite materials, has revolutionized 
structural design, bringing with it many advantages. Composites 
consist essentially of high strength -continuous or short- 
fibers dispersed in a matrix of relatively lower strength. They 
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can be tailored to suit a given application. This can be done 
by choosing suitable combinations of fiber/matrix materials and 
orienting the high strength, low density fibers along preferred 
directions. Features like reduced weight, higher strength and 
stiffness, directional properties and so on can be easily 
achieved through composites. 

However, composites are generally anisotropic in nature 
and hence introduce complicating effects like in-plane-shear 
coupling, bending-stretching coupling and so on. Nevertheless, 
the effects are partially offset by exploiting symmetry in the 
geometry of the structure. 

Laminated composite plate elements form an important part 
of any aerospace structure. A thorough dynamic analysis would 
reveal important vibrational and damping characteristics of 
such elements and this could form a vital input to the 
designer. This is necessary since the stresses induced due to 
dynamic loads can be an order of magnitude higher than static 
loads. 

A detailed review of previous research work on vibration 
analysis of plates reveals that most of the researchers' have 
used trigonometric or hyperbolic admissible functions i 
involving tedious integrations for the stiffness and mass 
matrix elements > in the Ray le igh-Ri tz method or have used the 
Finite Element method. Very few have attempted to use simple 
algebraic orthonormal polynomials. A much lesser number have 
applied the Ritz technique with simple algebraic admissible 
functions to the vibration analysis of Laminated composite 
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plates. 

In the present work, rectangular Laminated composite 
plates have been analysed for free vibrational characteristics. 
The assumed displacement field corresponds to that of the 
Classical Plate Theory. The Rayleigh-Ritz technique has been 
employed and the orthonormal bases for the Ritz series 
expansion of the displacements have been generated by the 
Gram-Schmidt procedure. 

[1.21 Literature Survey; 

In accordance with the conclusions drawn from the survey 
of available literature on plate vibrations, as has just been 
mentioned, it is only appropriate that the vast spectrum of 
research work in this field be divided into two distinct 
categories. They are as follows: 

¥■ 13 Dynamic analysis whSrein conventional admissible functions 
such as trigonometric and hyperbolic beam characteristic 
functions have been used in a Ritz technique or the Finite 
Element method has been employed. 

The work of Reissner and Stavsky^ provided the theoretical 
base for the study of vibration and buckling of laminated 
composite plates. They were one of the first few to consider 
the bending-stretching coupling effects in composites, in their 
analysis . 

2 

Bert and Mayberry used trigonometric and hyperbolic 
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functions for free vibration analysis of un-symmatr i cal ly 
laminated composite plates with all edges clamped. They have 
also presented corresponding experimental values for the first 
few frequencies. 

An early and notable effort in anisotropic plate analysis, 
galvinised by the coming of Composite materials, was that due 
to Ashton and Maddoups^. They demonstrated the use of beam 
eigenfunctions together with the Ritz method in vibration 
analyses of rectangular anisotropic plates. 

Whitney and Leissa presented closed form, exact solutions 

of free vibration frequencies and mode shapes for 

ant i -symmetr i c cross-ply and angle-ply laminates with 

simply-supported boundary conditions. Considering the 

classical types of boundary conditions (Clamped, 

Simply-supported or Free) that could exist on each edge of an 

un-symmetrical ly laminated plate, the number of independent 

problems possible, are many. In fact 16 different possibilities 

can be identified for each edge. One such possibility of a 

plate with simple supports <N .= 0 instead of u.= 8 in the 

nt t 

simply-supported boundary condition) was analysed by Whitney 

5 

and Leissa . They used a generalised Fourier series. 

6 

Later Whitney showed, using clamped boundary conditions 
for square anisotropic plates, how a generalised Fourier series 
solution led to an infinite frequency determinant different 
from that of the Ritz approach yielding more rapidly converging 
results. 

In an exce 1 lent paper on buckl ing and vibration , Jones^ 
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presented an exact solution and results for simply-supported 
plates, un-symmetr i cal ly laminated about their middle surface. 
Jones showed that bending-extension coupling due to asymmetry 
reduced vibration frequencies for common Composites like 
Boron/Epoxy and Graphite/Epoxy. It was also shown that for 
ant i -symmetr i c laminates, the effect of coupling died out 
rapidly with increasing number of layers but that for generally 
unsymmetric laminates the decay was slower .Simple trigonometric 
functions were used for the displacements. The largest effects 

were shown to occur for high modulus ratio materials. 

8 

Lin and King showed that exact solutions could be obtained 
for a class of unsymmetr ical ly laminated plates having two 
opposite sides simply-supported with arbitrary edge conditions 
on the other two sides, by suitably choosing the admissible 
functions. 

9 

Bert derived a relation to compute frequencies of a 
generally anisotropic plate having arbitrary boundary 
conditions, given the frequencies of a corresponding 
rectangular isotropic plate. Later, the same investigator used 
this result and suggested an optimal design procedure to 

maximise the fundamental frequencies of such plates. 

10 

The work by Bert and Ghen, aimed at studying the effect 
of shear deformation on vibration of anti -symmetric angle-ply 
laminates. They obtained closed form solutions for all edges 
siit^ly-supported, using the YNS < Yang-Norri s-Stavsky > theory. 
They too, employed trigonometric functions for the displacement 
and slope variables. A high modulus (Graphite/Epoxy> mater 
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with Ej^/E^ = 40 was used and results were presented showing the 

effects of plate aspect ratio. length/thickness ratio, 

lamination angle and number of layers. It was also shown that 

neglecting in-plane and rotatory inertias had little effect on 

the fundamental frequency, even for thick plates. 

Dickinson^^ initiated the use of the Simply-supported 

plate functions which were described as admissible functions, 

derived from the mode shapes of vibration, of plates having two 

parallel edges simply-supported, the boundary conditions on the 

other two edges being chosen appropriately for the plate under 

consideration. It was shown that accuracies resulting from the 

use of these functions were superior to those from the use of 

beam functions and the more complicated Bolotin functions. 

1 2 

Reddy, used the YNS theory in his paper on vibration 

analysis of anti -symmetric laminated plates which was a close 

1 % 

parallel to the one published by Bert and Chen .However, here 

the Finite Element method was employed. The plate was 

discretized using eight noded Serendipity elements. Results for 

30° and 45° angle ply laminates for two different modulus 

ratios’ (typical of Graphite/Epoxy) with Ej_/E.j. = 40 and 25 were 

1 9 

pr€isiifit#d and cofnpal^ed with those obtained by Bert and Chen 
Also frequencies obtained for a simply-supported thick 
Isotropic plate were compared with those due to 5-D linear 
elasticity theory and that due to Hindi in’s. 

In a later paper, Dickinson and Li^^ used the previously 
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described simply-supported plate functions to analyse flexural 
vibrations of rectangular plates by the Ray le igh-Ri tz method 
and confirmed that improvement in accuracy of the natural 
frequencies was indeed possible, compared to beam functions. 
However. it was also concluded that 'Simply-supported' 
functions were less satisfactory for plates with free edges. 

Bhimaraddi and Stevens^^ proposed a two dimensional theory 
for the analysis of thick homogeneous laminated plates. The 
cubic order polynomial form for in-plane displacement 
quantities was maintained ensuring at the same time the more 
realistic parabolic variation for transverse shear strains. 
Further, trigonometric displacement functions were used in the 
analysis. 

In another lucid and excellent paper, Lin et al^^ 

presented their research work with the objective of predicting 

the natural frequency and Specific Damping Capacity <SDC> of 

laminated Composite plates in various modes of vibration by the 

Finite Element method. An 8-noded, 40 degree of freedom element 

was used. The damped element model was developed by introducing 

^1 '^2 beam specimen tested in flexure- 

and -the SDC of 0*^ and 90° specimen tested in 

12 25 

longitudnal shear- in the strain energy expressions. 

Experimental results were also presented and the 
measurement technique was based on the transient testing 
technique using a zoom-FFT. A micro-computer interfaced with 
the FFT analyser recorded and calculated the frequency spectrum 
using resonant curve fitting or half power point technique to 
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obtain the natural frequency and vibration damping immediately. 

Malhotra et al studied the effects of fiber orientation 
and boundary conditions on the frequencies of thin square 
orthotropic plates. Comprehensive results for various boundary 
conditions and lamination angles were presented for four types 
of Composites: Glass, Kevlar. Boron, Graphite/Epoxy. 

► Recent work on study of free vibrations of plates using 

simple algebraic orthonormal polynomials; 

1 7 

Bhat made use of algebraic functions in plate vibration 
studies. Bhat obtained the natural frequencies of rectangular 
isotropic plates using a set of beam characteristic orthogonal 
algebraic polynomials in the Ray le igh-Ri tz method. The 
orthogonal polynomials were generated using the Gram-Schmidt 
procedure. The first term was so chosen as to satisfy all the 
boundary conditions of the corresponding beam problems which 
were a subset of the plate boundary conditions and the 
succeeding terms satisfied the geometric boundary condition 
only. Thus the basis for the series expansion of the 
displacements were formed of products of the x and y polynomial 
terms so generated, with the procedure automatically ensuring 
completeness and continuity of these polynomials. 

Excellent results for plates with simply-supported, 
clamped and free edge conditions were presented which compared 
very well with the exact value and with the results obtained by 
other investigators. 

1'8 ■ " ' 

Dickinson and Blasio made use of the orthogonal 
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polynomial functions proposed by Bhat^^ for vibration and 
buckling analysis of isotropic and orthotropic plates, 
employing the Ritz technique. 

The authors confirmed that these simpler polynomials 
provided an attractive alternative for admissible functions as 
compared to the conventional trigonometric and hyperbolic 
functions, such as those used in the ’Simply-supported plate 
functions’ or the ’Degenerated beam functions'. 

It was proposed that in the case of plates with free 
edges, simpler starting terms (of lower orders) could be used 
with very little loss in accuracy. Vibration and buckling 
problems of isotropic and specially orthotropic plates, 
subjected to in-plane direct and shear loading were addressed. 
Finally, it was concluded that the simple polynomials relaxed 
the over restraint in the case of plates with free edges and 
rendered the product integrals for the stiffness and mass 
matrices, easy. 

IQ 

Bhat proposed a further modification to the approach 
using algebraic polynomials. Orthogonal polynomials in two 
variables were used with a little difference in the earlier 
approach (Bhat^^> and the new one. These polynomials were used 
to obtain the natural frequencies of triangular plates of 
different configurations. 

Lam et al^® presented results of vibration analysis of 
single layered composite plates. They used a modified version 
of Bhat’s two dimensional orthogonal plate functions, in a Ritz 
analysis. The effects of boundary conditions, aspect ratio and 
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fiber orientation on the natural frequencies of Composite 

plates Mere studied. 

21 

Liew et al used orthogonal plate functions for vibration 

analysis of rectangular isotropic plates including those with 

free edges. A starting term similar to that of Bhat^^ was used 

but with a modified Gram-Schmidt generation of succeeding terms 

which led to faster convergence. Very good results were 

obtained for plates with various types of boundary conditions. 

22 

Liew and Lam used the orthogonal plate functions (two 
dimensional) in a Rayleigh-Ritz approach to investigate free 
vibration of isotropic and anisotropic trapezoidal plates, with 
simply-supported, clamped and many other types of boundary 
conditions. 

[1.31 Objective and Scope of the Present Work: 

The literature review clearly indicates that the technique 
of using orthogonal algebraic polynomials is an elegant and a 
very powerful one for the vibration analysis of plates. 
Further, it is also evident that this technique has, to date, 
been applied only for the analysis of isotropic and orthotropic 
plates. Keeping this in view, an application of this technique 
to the dynamic analysis of general Composite laminates has been 
attempted in the present work. 

The present work focuses on the free vibration analysis of 
Composite laminates with different boundary conditions. 



11 


materials, fiber orientation and lay-ups. The Classical Plate 
Theory has been used and orthogonal polynomials have been 
employed in the Rayleigh-Ri tz analysis. A computer program <in 
Fortran 77> that generates the orthonormal polynomials 
automatically and computes the necessary product integrals has 
been developed. Simple Gauss Quadrature and numerical 
differentiation schemes have been used. Eigenvalues and 
eigenvectors have been computed using available HAG routines 
all these on the HP 9000/lB50s . 

Results for two types of materials have been generated for 
antisymmetric angle-ply and cross-ply and symmetric cross-ply 
laminates for various boundary conditions. The usual check for 
the isotropic case has also been made. The results have 
compared with those generated by standard techniques such as 
Finite Element techniques and so on. Some results for clamped 
isotropic tapered plates have also been included. 

Chapter 2 covers a detailed description of the problem 
formulation. 

Chapter 5 contains an analysis of the results presented. 

Chapter 4 contains conclusions drawn and suggestions for 


future work. 
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CHAPTER 2 

THEORY AND FORMULATION 

12.1] General Solid Continuum analysis: 

[2.1.11 Degrees of freedom of a system: 

The minimum number of independent coordinates required to 
completely determine the configuration of a system. 

[2.1.2] Material Point; 

A point in the domain of the continuum where a mass point 
is assumed to be located. 

[2.1.31 The Solid Continuum; 

A solid continuum is assumed to be a continuous material 
that is devoid of voids and gaps. The configuration of a 
deformable solid continuum is determined by the simultaneous 
positions of its infinite material points relative to a fixed 
coordinate system. Thus a deformable solid requires 
specification of an infinite set of degrees of freedom. each 
corresponding to the position of a material point. 



1 > 

Ideally, the position of a material point is uniquely 
determined by three independent coordinates. However, it is not 
practicable to deal with discrete points and hence 
infinitesimal volumes are considered. These are characterised 
by one or more field variables such as displacements, slopes 
and rotations. In turn, each of the field variable is expressed 
as the sum of an infinite series of linearly independent 
functions. For example, if A is a field variable, we have 

A = Z c. £.(x,y,z) <2.1> 
i=l ‘ ‘ 

The linearly independent basis set spans a 
configuration space of the system and any subset of it spans a 
subspace. 

In a Dynamic analysis, such as in the free vibration 
problem, the contemporaneous positions of the material points 
relative to a fixed coordinate system describe its 
configuration. 

12.1.41 The Problem; 

The objective of any analysis using various techniques 
such as those described below is to solve for the unknown 
coefficients c^ , from equations of the type 
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[K] (c) = (f 3 

for a static problem, or 

[K] CO [M] CO 

for a dynamic eigenvalue problem, where 

|[ K 3 is the stiffness n^trix 

[ M ] is the mass matrix 

[ f 3 is the load vector, and 
X the eigenvalue 


The coefficients c^ can be visualised as components of the 
solution vector in the configuration space for the given set of 
loads on the system. Thus they determine the extent to which 
each of the degrees of freedom, participate in the 

solution. 

ZZ.1.31 Some Analytical techniques: 

It is obviously impossible to solve for an infinite set of 
coefficients and hence many of the techniques represent the 
field variables as the sum of a finite number of terms <say N> 
of a series as, 

N 

Z c . s . 

I 1 


A 


( 2 . 2 ) 



► a) The Rayleigh-Ritz technique: 


In the Ritz analysis, each field variable of the problem 
is expressed, as the sum of a finite number of terms of a 
series. 

The total potential energy, V and the kinetic energy, T, 

of the system and hence the Lagrangian function (L = T - V) , 

computed by using equations of the type <2.2> in the respective 

energy expressions, are rendered functions solely of the 

unknown coefficients, c.. Thus we have, 

1 


V 

T 

L 


V(c, , c_ , . . , c.,) 

T<Cj .c^, . • ,Cj^> 
L-CCj *^2* * * * ^ 


(2.5> 


Invoking the Principle of Stationary Value of the total 
Potential and taking a variation on V, we have for a system in 
equilibrium. 


Xt, f , X . 

<5V = I — - 6c. + — ^c^ + 
^ dCj dc^ 


. . + — o c. 


9c 


M 


h) 


( 2 . 4 > 


Now since the 6cj are all arbitrary and independent, they 
are all not necessarily equal to zero and so we have. 
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— = 0 . i = 1 .2. . . .M (2.5) 

9c 

1 


Similarly for a system in Dynamic equilibrium undergoing 
Harmonic oscillations, we can write 


9 S 

t. 


L dt 


9 c . 


= 0 


< 2 . 6 ) 


To ensure that the Ritz approximations result in a 

solution and that it converges to the true solution as the 

number of terms of the series are increased, the choice of the 

set C£^ . > has to satisfy certain conditions are detailed below: 

1 

# i) Completeness: The choice of the s functions should begin 
from the lowest order of the polynomial being used, preferably 
from the zeroth order and can increase upto any desired order 
without intermediate orders being missed. 

The effect of ignoring certain intermediate orders, 
provided the first few terms are complete, is to bias the 
solution towards certain preferred modes. 

Completeness can easily be ensured by choosing complete 
sets of polynomials from the Pascal’s triangle shown below. 
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• ii> Boundary Conditions: All the e. chosen should satisfy the 
Essential boundary conditions of the problem. They may or may 
not satisfy the Natural boundary conditions. 

• iii) Continuity; The functions chosen should be continuous at 
least upto the order required by the variational statement 
being used ( usually C^ continuous). 

The £ ^ fulfilling the above requirements constitute an 
"Admissible Set”. The Ritz technique has been used in the 
present work. 

► b> The Galerkin technique: 

In this technique, too, the variables are expressed in the 
form given by eq<2.Z>. However, here, the Differential equation 
is used instead of the variational statement. 

If £. denotes an operator, we have 

M 

= L< S c. s.y 
i=l * 


R 


( 2 . 8 ) 
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where R is the error or residual due to the finite 
approximation. This error is now minimised by forcing it to be 
orthogonal to each of the basis thereby reducing each of 
its components along to a minimum. 


i.e jRfi.dv = e 1=1, 2, ..,M ..<2.9> 

V 

The completeness condition as mentioned previously in the 
Ritz analysis applies here also, in choosing an ’Admissible 
Set'. However, the functions have to satisfy all the boundary 
conditions, essential and natural. Continuity requirement is 
also generally higher and has to correspond at least to the 
order of the highest derivative of the Differential equation 
used. 

Normally the choice of the basis set Ce ^ } is arbitrary. 
However, an Admissible Set of orthonormal functions can be very 
elegantly generated using the Gram-Schmidt procedure. This also 
ensures all the requirements of completeness, boundary 
conditions and continuity automatically. 

[2.1.61 The Gram-Schmidt Procedure: 

Given a set of linearly independent vectors; 
lj(s>,l2<s>,..., l^(s> 

one can construct a set of orthonormal vectors 
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9j<s> .g 2 <s) 

by the Gram-Schmidt procedure as folloMs 

1 


•l 


1 

IMill 


= 1 


<92 - Oj 


l|p,ll 


Pj = Ij - <gj.lj> gj - <g 2 . 1 j> g^ 


= 



IIP3II 


and so on, where 


I Ml I I ~ norm of Ij 


<g 2 -l 2 > = scalar or inner product of vectors gj and 1^ 

^2 




[2.1.6.11 A Geometric interpretation; 

For M <= 3, an interesting geometric interpretation can be 
easily given to the Gram-Schmidt procedure. AM the vectors 
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reduce to the common three dimensional vectors of the 
Euclidean space. 

Refering to fig <2. 2a), by subtracting from l^. its 
projection on g^ , we obtain the orthonormal vector g^. Vector 
g^ is now normal to g^ . Here. 

I Ml I I gives the length of vector 1^, and 

is the cosine of the angle between g^ and l^* 

A similar argument can be extended to explain fig (2.2b> 
for three dimensions. 


12.2] Laminated Composite Plate Analysis : 

A plate is a solid continuum with a regular geometry and 
with one of its dimensions being very small compared to the 
other two lateral dimensions. 

(2.2.1] Classical Plate Theory; 

The Classical Plate Theory < henceforth called CPT > has 
been used in the present analysis. Plates are conventionally 
classified as thick plates <a/h < 20) or thin plates <a/h >= 
20). CPT is the small deflection theory of bending of thin 
plates and it involves the following assumptions. 

► The mid-surface displacements are small compared to the 
thickness of the plate. 
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► The mid-surface of the plate is unstrained and is the 
neutral plane after loading. 

► Plane sections normal to the mid-surface remain plane and 

perpendicular to the midplane ( implies negligible 

transverse shear strains e and e >. 

x 2 yz 

► o can be neglected. 

z 

In accordance with CPT, the following displacement field 
has been assumed. 

dvi 

u<x,y,z> = u <x,y> - z — 
o dx 

dw 

v<x.y.z> = v^<x.y> - ^ ••• <2. 11) 

w ( X , y , z > = w^ < X , y ) 


(2.2.21 Strain Dispalcement Relations: 


Using 

the 

above field. 

we have 

the 

following 

relations for 

strains 







e = 

X 

^u 

dx 

du 

o 

dx 

2 

d w 

o 

^ 2 
ax^ 




It 

dv 

dy 

1 

O 

> 

„2 
a w 

o 

2 

^y 



<2.12> 


du 

a au 

av 

a w 

o 


^xy ■ 

3y 

^ ax ay 

ax ^ 

dx 

9y 




zz 


[2.2.3] Stress-Strain Equations; 

A laminate is made up a number of laminae. Each lamina 
has, in general, fibers oriented at a given angle. Thus the 
individual laminae are orthotropic in their own plane with the 
principal planes of orthotropy being oriented along and normal 
to the fiber direction. 

If L and T stand for the fiber direction and its normal, 
we have the following constitutive relations for the 
lamina , 



< 2 . 1 >> 


Since the L/T directions vary with each lamina, and that 
in general they do not coincide with the structural axes, the 
stresses and strains of each lamina will have to be transformed 
along the structural axes. Following the procedure given in 
Ref. tZ61 the transformed constitutive equations in terms of 
the reference axes is given by 
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[2.2.4] Stress Resultants: 

The laminate stress and moment resultants, as shown in 
Fig. <2. 5), are defined as follows: 

h/2 


N N N 
X. y, xy 


1 


ia ty T ) dz 
X. y. xy 


-h/2 

h/2 


M M M 
X, y, xy 


I 


(O' o T ) Z dz 
X , y , xy 


-h/2 


(2.15) 


Substituting relations <2.14> in <2. 15). we have 




N 


N 


xy 


M 


M 


} - 


’ll 


^12 ^22 


^16 ^26 


sym 




11 

®12 

®16 

°11 


12 

®22 

®26 

°12 

°22 

16 

®26 

®66 

°16 

°26 


-1 

X 

o 
e 

y 

o 

y' 

xy 


. 


n 




16 °26 *^6^ 


<2.15a> 


where, the x's are curvatures and. 

2 . 


NL p 
(A . B. . D. .) = E 

n=l 2 


■n+1 


Q^"^C1. Z. Z^)dz 


(2.16) 


n 


NL : Number of layers in laminate. 
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[2.2.51 Potential and Kinetic energies of the laminate: 


In the absence 

of 

applied loads 

and 

body 

forces , 

the 

potential energy of 

the 

plate, in terms 

of 

the 

stresses 

and 


strains, can be written as follows: 


a b +h/2 




0 0 -h/2 


(o' e + O' B + O’ e 
X X y y z z 

+T y +T r +T y ldxdydz......(2.17) 

xy xy xz xz yz yz-' 


where 

a = length of plate 
b = width of plate 
h = thickness of plate 

Non-dimensional ising the variables x and y with respect to 
the plate dimensions, we have 


? = x/a 

17 = y/b 

Using eqns. (2.14) to <2.15a> and <2.17), we obtain the 
potential energy in the non-dimensional form as below: 
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du 


du du 


2a ^'Aj 2 [ j + 2a j + 2a Aj^|^ j + 


du dv 


dv dv 




5u 

av .. 

0 

° 1 

dK 

W J 

. au 

a^w 

0 

1 0 


2 

a?^ 


2 

^v dw 


du d w 


-T r ^ ^ VW -v *% rf- W4 V w -^. 

2 a M B 22 [ ^ ^ ] - ^ j 


2 2 
-dvdw^ - f d\i S Vi ^ 

®66 [ 5r 5f4) - ®i2 [ -;;;|- ] 


By B ^ 

®12 ( TT J ■ ®16 


^ du dw Vi ^ 

fz — ° °- + — ° ° 1 

d? a? dT7 ^ dr? ^^2 J 


2a B 


. av a w - , 

16 ( &r ] - 2a^ B 


au_ a w 


a^ 




2 2 

- av a w av aw 


-2 
^ a w 

..2 . 

2 

r ^ w 

r, 0 0 0 o') 

1 + n 

f 

J ^ °22 

r ° 1 

ar? a? ar? a? ^2 J 

1 ^ *^11 

1 a?2 

i- ar?" J 


a Vi . 2 


«2 ^2 
aw aw 


^ ^ °66 ( a? ^ ] ^ ^ '^12 ( ■”^;^ ) 


a2 -2 

( aw aw 

— 5 ^ ^ 

£ CF< J ^ ^ 


.2 ^2 

d w B ^ 


V V ■» J- V W V 

-^ ] + 4 /J 0^6 [ ] 


d? dr? 


C2.18> 


d? dr? : differential area of plate mid-surface 
?j = a/b = plate aspect ratio. 
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The kinetic energy of the plate in its non-dimensional 
form is obtained as; 


1 1 +h/2 


pab 

2 


HI 

9 9 -h/2 


[ <u>^ + <v>^ + <w>^ ] d? d77 dz 


<2. 19> 


where, p is the density of the material. 


The Ritz series expansion of the field variables in eq. 
(2.11) was carried out with the following orthonormal 
functions: 


.T?> 

.ri> 




+ C 


u 

12 




+ 





+ C 


12 




f" .... 





+ C 


w 

12 




+ . . . . 


< 2 . 20 ) 


The following recurrence formulae as proposed by Bhat 
which are based on the Gram-Schmidt procedure, were used to 
generate the above functions: 

= (? - b_> a- 

? 2 • ^ 1 

a- = (? - b,) 

^f ^ ^f-1 


< 2 . 21 ) 




f-2 
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In general, a pair of constants <b^, c^> are required to 
build the f term and they are determined using the following 
orthogonality condition 

- 0. if i j 
^1. if i = j 



which finally leads to the expressions for the constants as 
given below: 


b 


f 






2 

«« cl? 

^ f-1 


oi^ = constant x a„ 

? f Kf 


where the constant is a normalising factor and is given by 


constant 


1 



d? 


The same set of equations (2.21) were used to generate the 
?,r) functions for all the variables. These functions depend 
solely on the boundary conditions being used. 
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[2.2.61 Boundary Conditions: 

The following boundary conditions were used in the present 
anal ys i s : 

• i) For Anti -symmetric cross-ply laminates 

S2 ; w = 0, M =0. N =0, u. = 0 

n n t 

CZ ; w=0, w = 0, N =0, u. =0 

, n n t 

• ii> For Ani-symmetric Angle-ply laminates 

S5 ; w = 0, M = 0, u = 0. N . = 0 

n n n’t 

C5 : W = 0. W „ = 0. = 0. N^. = 0 

, n n nt 

where. 

S => simply - supported 

C => clamped 

n => direction normal to edge being considered 

t => direction tangential to edge 

,n => derivative with respect to variable along direction n. 

One among the many advangates of this method is the 

relative case and elegance with which starting terms for the 

Ritz series may be chosen. The starting terms of all variables 

satisfy both the geometric and the natural boundary conditions 
whereas the succeeding terms satisfy only the geometric 
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boundary conditions: 


The following starting terms were chosen for the boundary 
conditions indicated: 

S2. C2 : 

Z 

a„ = constant <f - F ) 

(3 = constant <0.006 - 


ftj, = constant <0.006 - 2?^ + 5?^) <2.22> 

^ 1 

2 

a = constant <77 - t) > 

S5, C5 : 

a and ft were interchanged to give S5, C5 boundary 
conditions. 


The 4> functions were chosen on the lines of Bhat^^. These 
functions satisfy the equivalent beam boundary conditions <both 
geometric and natural). 


S : 4>j, = constant <? - 2?^^ + ?^) 

^ 1 

Z I 4>f - constant <?^ - Zt} + 

^ 1 


<2. 23) 


The functions are identical to the functions with 77 
replacing f. All the constants above are the normalising 
factors of the respective terms as explained earlier. 
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12.3] Hamilton's Principle: 


Hamilton's Principle is a generalisation of the principle 
of virtual displacements to the dynamics of systems of 
particles, rigid bodies and deformable bodies. The 
Eul er-Lagrange equations derived from this principle result in 
the familiar Newtonian equations of motion for a generalised 
set of co-ordinates. 

The principle states that 'of all possible paths a system 
could travel, from its position at time tj to another position 
at time t^ , the actual path is that path which renders the 
integral 


J 


L dt 


an extremum. That is 


61^ = <5 J L <qj, q^. qj, q^...) dt = 0 

t , 


<2. 24) 


where L is the Lagrangian function. From eq <2. 24) we also get 
the Eul er-Lagrange equations as follows: 


ah 



= 0, i=1.2,5.. 


<2.25) 


q^^ are the generalised coordinates 
dot represents derivative with respect to time. 
For synchronous motion due to free vibration, we have: 
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u = 

V = 

w = 


u e 


V e 


M e 


jwt 

jwt 

jcot 


(2.Z6> 


where the system vibrates with a natural frequency to. 

Using equations (2. 11). (2.18), <2. 19), <2. 20) and (2.26), 
we obtain the Lagrangian function L, as 

L = T - V 


Taking a variation on the set of unknown coefficients we have. 



L dt 


0 


<2. 27) 


Factoring the time integral f e dt out of the 

above equation and performing the necessary differentiations 

with respect to the unknown coefficients, we arrive at the 
familiar eigenvalue problem 

[ K ] { C f = X [ M ] -{ C } ... <2.28) 

where X is the eigenvalue and the non-dimensional frequency 
parameter is given by 

i"X = CO a^ ... (2.29) 

K and M are the stiffness and mass matrices respectively and C 
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stands for the vector of unknown coefficients as also the eigen 
vector . 

The number of terms (denoted by 'noft'), used in the 
series of eq (2.20> for all the variables, in both the f and rj 
directions are the same. Thus the final size of the stiffness 
and mass matrices can be shown to be 5 x noft^ . 


(2.3.1) Components of the Stiffness and Mass Matrices; 


The following type of notation used for the various 
sub-matrices that go to build the stiffness and mass matrices 
needs some explanation. 

J . [ »* p* ] . [ r v-" ] 

choosing a value of noft = 2 (i.e. two terms) to explain the 
above notations, we arrive at the following 4x4 matrices 




n 

.m 

n j 

^1 

^1 


-do- 


n 

.m 

n 1 

^2 


^2 

'"2 1 

- 

— do 


T 

1 -do- 


d? dr? 


All sub-matrices occur as products of a pair of matrices 
as shown above. The dot indicating the product DOES NOT stand 
for the usual matrix multiplication but implies that the 
corresponding individual elements of the two matrics should be 
multiplied. Further, all the elements of the first matrix are 
always, functions of the variable ? only and those of the 
second matrix that of t? only. The super scr ipts i , j , m and n 
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stand for the order of differentiation of that particular term. 

2 2 

All the submatrices are of size noft x noft and the 
sub-submatrices are of size noft x noft. 

The stiffness and mass matrices can now be expressed in 
terms of the various sub-matrices using the notation explained 
above. So we have for the stiffness matrix. 



■[SjjJ 

[S,^] 

[Si,] ■ 

[ K] = 


[ 

[S„] 


SYM 


[S„] . 


where 

[Sjj] =CAj,[al] +Aj^[a6J tAjJal6J> 

[Si^J "CAj^taU] +A^^[a6^^] tAj^[al6^^] + A^^t aZ6^^ ) 

[ S^J,] = CAjjC aZ) + 

[Sjj] =(Bjj[bl] ^Bj^fbUJ t Bj^[ bl6J + J b26J ) 

[Sj,] "CBjjCbZJ +BjJb6^] +Bj2[bl2^] t Bj^[ bl6^] + B^^L b26^] > 

[Sjj] =COjjtdl} +D^^lciZ} +D^Jd6] + 

+DjJdl63 +D2^[d26]> 

Details of all of the above matrices are given in the 
Appendix. 

The mass matrix may be expressed as follows; 
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[M] 


where [ I] 

[ Mj] = [ i; 


[ I] 

8 

8 


[ 1 ] 

8 


SYM 


[Mj] 


stands for the identity matrix, and 

I ] 






FIG. 2-1 CONCEPTUAL REPRESENTATION OF SOLUTION VECTOR 
's' IN THE CONFIGURATION SPACE 




(°) 



FIG. 2-2 GEOMETRIC INTERPRETATION OF THE 
GRAM- SCHMIDT PROCEDURE IN 
2-D AND 3-D SPACE 





Fig.2.3 STRESS RESULTANTS OF A 

LAMINATE 
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CHAPTER 3 


RESULTS AND DISCUSSIONS 


[3.1J Introduction; 

The governing equations derived in Chapter 2 have been 
solved for free vibration analysis of composite plates. A 
computer program has been developed on the Mainframe HP 
9000-<^50s using available NAG routines for the eigenvalue 
computations. 

The core of any dynamic analysis < such as the present 
free vibration analysis > lies in the generation of the various 
stiffness and mass matrices. However, in the present analysis, 
automatic generation of the Ritz series expansion functions 
formed an equally important part, being inputs to the stiffness 
and mass elements. 

A computer program incorporating the function generation 
aspect was developed to analyse out of plane vibrations of an 
isotropic plate and the results were checked with those 
available in the literature. The program was then modified to 
suit product integrals involving cross products of functions of 
different displacement variables. Sub-routines to generate the 
A.B and D matrices for the composite plate were also developed 
to produce the final program. This program was rechecked for 
numerical accuracies when reduced to the isotropic case and was 
then used for composite plates. The boundary conditions were 



59 


incorporated by changing only the starting terms of the series 
to obtain results for the new conditions. Finally the program 
was further modified to obtain results for an isotropic clamped 
tapered plate. 

In order to perform the necessary integrations and 

differentiations for generating the various matrix elements. 

Gauss Quadrature and Numerical differentiation schemes were 

employed. The number of Gauss points used in the Gauss 

Quadrature integration scheme determine the highest order of 

the polynomial that can be integrated accurately. For example. 

if ng is the number of Gauss points used, a polynomial of order 

< 2*n-l can be integrated very accurately. In the present 

analysis, the order of the polynomial and the number of terms 

are related, i.e., with every new term, the increase in the 

order of the polynomial is by one. Consequently, the number of 

terms and the number of Gauss points also get related. Only 4 

terms were used < number of terms = noft = 4 > in the analysis 

as it was found that this was sufficient for convergence of the 

first three to four frequencies. Ten Gauss points were used for 

the integration. The total number of integrations for the most 

general laminate ( i.e., assuming that the A.B and D nraitrices 

2 

are all fully populated > would be 70 * noft . The average 

total time for execution of one run of the program was less 


than a minute. 
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13. 2J Results: 


The composite materials used in the present analysis are 
typical of High and Low modulus Graphite/Epoxy variety. Their 
properties have been presented in Table 5.1. 

[3.2.11 Study of Convergence: 


The size of the stiffness and mass matrices is given by 
2 

3*noft . As mentioned previously, all results have been 
obtained using noft = 4, This leads to a final matrix size of 
48 X 48 for the stiffness and mass matrices. 

Convergence studies of the non-dimensional frequency Jx 
were carried out for isotropic, orthotropic and composite 
plates. Figure 5.1 shows a typical behaviour for the case of a 
simply-supported, square, anti -symmetric angle-ply laminate. It 
is evident from the figure that 4 terms are sufficient to 

^ I*) 

ensure convergence upto the 4 mode. It was also observed that 
the rate of convergence exhibited slight oscillations. 

[3.2.21 Effect of Side/thickness ratio: 

Figure 5.2 shows the variation of the non-dimensional 
frequency with various side/thickness ratios for a 
simply-supported, square, ant i -symmetric angle-plY laminate. 
The near flat response seen, is to be expected, since the 
Classical Plate Theory neglects transverse shear effects. It 
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generally over predicts the response but is fairly accurate for 
thin plates < a/h >= 20 >. Thus for thick plates, -where shear 
effects are predominant- the results are conservative. 

[3.2.31 Isotropic and Orthotropic plates; 

Table 5.2 shows a comparison of the non-dimensional 

frequency for a simply-supported, square isotropic plate. 

Solving for these out of plane frequencies required a value of 

noft = 4, for convergence of the first 4 frequencies but the 

matrix sizes were only 16 x 16. As can be seen from the table, 

the results by the present method are in excellent agreement 

17 21 

with those of Bhat and Liew . Frequencies for various aspect 

ratios were also checked. It n^y be pointed out that the the 

choice of the starting functions and the generation of the 

succeeding terms is much simpler in the present method as 

21 

compared to that due to Liew . The constant 0.006 in the 
starting terms of the Ritz series shown in Chapter 2 needs 
explanation. The constant was introduced to prevent the 
complementary in-plane displacement from going to zero at the 
specified edge. The numerical value of 0.006 for the constant 
was arrived at after a number of trials with various other 
values were carried out, for minimising the frequencies. 

In Table 5.3 are shown frequencies of an orthotropic, 
square, simply-supported plate obtained by employing the same 
functions that were used for the isotropic plate studies. It is 
observed that the results obtained by this technique are in 
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close agreement with those seen in Ref [253. 

( 3« 43 Symmetric Cross-ply laminates: 

Figures 5.5 and 5.4 show the variation of the 
non-dimensional frequency, -jx with aspect ratio for four 
layered symmetric cross-ply laminates with simply-supported and 
clamped boundary conditions along the four edges, respectively. 
The frequencies for the clamped case are seen to be higher than 
the corresponding frequencies for the simply-supported case. 
This is due to the fact that the clamped boundary conditions 
tend to increase the overall stiffness of the system thereby 
requiring higher energies to set it vibrating. The frequencies 
are also seen to be increasing with aspect ratio. This 
indicates that the geometry of the plate also has an influence 
on the frequency. Physically it maybe explained as follows. The 
fundamental or lowest frequencies of a plate can be attributed 
mainly to the out of plane vibration modes. This is because the 
in-plane stiffnesses are higher compared to the flexural 
stiffness and hence the in-plane frequencies also higher. Now, 
the sub-matrix for out of plane stiffness is a direct function 
of the aspect ratio. Hence it is evident that for low aspect 
ratios', this dominant sub-matrix < dominant as far as the 
lower frequencies go > has a fairly low value and the 
corresponding fundamental frequencies are also low. Naturally 
the frequencies would increase for higher aspect ratios. It is 
also seen that the frequencies for the low modulus material are 
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lower compared to that of the high modulus material. This is 
simply because the plate of the higher modulus material is 
st i f f er . 


(3.2.51 Ant i -symmetr i c Cross-ply laminates; 

The variation of non-dimensional frequency with aspect 
ratio for ant i -symmetr i c cross-ply laminates with 
simply-supported and clamped boundary conditions are shown in 
figures 5.5 and 5.6. The frequencies for this case are seen to 
be lower compared to the ones for the symmetric cross-ply 
laminates. This behaviour maybe explained by considering the 
fact that for symmetric laminates, the elements of the B matrix 
are all zeros. Since the B matrix couples in-plane and bending 
displacements, it tends to make the plate more flexible. Thus, 
when the B matrix is zero. the stiffness is increased and 
consequently frequencies also increase for such laminates. The 
variation with aspect ratio and material is seen to be similar 
to the previous case as the behaviour is the same as explained 
earlier. 

(3.2.61 Ant i -symmetr ic angle-ply laminates; 

Figures 5.7 to 5.12 show the variation of the 
non-dimensional frequency with aspect ratio for anti -symmetric 
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angle-ply laminates with fiber orientation of Ijtje^and 4?” and 
for simply-supported and clamped boundary conditions. The 
nature of variation of frequencies with aspect ratio and 
boundary conditions are seen to be the same as in the case of 
cross-ply laminates and can be explained in a similar manner. 
However it is interesting to study the frequencies of plates 
with aspect ratios < 1 having different fiber orientations. It 
is seen that for low aspect ratios the frequencies decrease 
with increase with fiber orientation angle i.e., . 
This is because for a given value of aspect ratio, there is a 
preferential increase in stiffness along the shorter length of 
the plate for small fiber orientation angles, which makes it 
more difficult to excite such plates. However for fiber 
orientation angle of 45*. no such preferred direction exists 
and the structure is more flexible. But for aspect ratios > 1 
the trend reverses. 

Figures 5.13 to 5.16 show the variation of non-dimensional 
frequency with fiber orientation angle for simply-supported and 
clamped plates of different materials. It is known that as the 
number of layer are increased the bending-extension coupling 
(due to the b,, and b,, elements of the B matrix > decreases. 
This increases the frequency of vibration. In figures 5.13 and 
3.14, it is clearly seen that for plates with more than two 
layers, the frequency increases with fiber orientation angle as 
expected. However, for two layered plates there is a decrease 
due to the dominating effect of coupling. In figures 3.15 and 
5.16, the decoupling effect is evident with the increase in 



number of layers, but the trends for the differently layered 
plates are all the same due to the clamped boundary condition. 

Table 5.4 shoMS a comparison of non-dimensional 

o 

frequencies for ant i -symmetr i c angle-ply laminates with 45 
fiber orientation. This table reflects the superior behaviour 
of the present method. It is seen that the results obtained are 
in excellent comparison with those of Ref [123. It would be 
appropriate to note that the results presented in Ref ZlZl were 
through an FEM analysis that included shear deformation. In 
both the cases the results have been obtained for a/h = 50, 

which is close to thin plates. The results indicate that the 
present method is quite sufficient for the study of slightly 
thick plates. 

13.2.71 Tapered plates: 

The isotropic analysis was modified to include a tapered 
plate with the taper on the thickness being along X direction 
only. A clamped boundary condition was used. Table 5.5 shows a 
comparison of frequencies for a square plate with different 
taper parameters. The taper parameter chosen is the slope of 
the taper. It is seen that the results by the present method 
match with those of RefC24] very well ( within 4* error >. The 
present method lends itself very well to a Ritz type of 
analysis with simpler starting terms as compared to the 
Galerkin. used in RefC243. wherein the starting terms are 
constrained to satisfy ail the boundary conditions. 

central 
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TABLE 3.1 


PROPERTIES OF COMPOSITE MATERIALS USED 


Properties Material-A Material-B 


E 


L 


E 


T 


48 


2 ^ 



E 


T 


0.5 


0.5 




LT 


0.25 


0.25 



TABLE 3.2 


COMPARISON OF NON-DIMENSIONAL FREQUENCIES 



FOR A SIMPLY-SUPPORTED SQUARE ISOTROPIC PLATE 


Mode 

Exact Value^^ 

17 

Bhat 

Liew et al^^ 

Present Method 

1 

19.759 

19.759 

19.74 

19.758 

2 

49.548 

49.548 

49.55 

49.547 

5 

49.548 

49.548 

- 

49.547 

4 

78.957 

78.957 

79.05 

78.956 

5 

98.696 

99.504 

99.25 

99.505 

6 

98.696 

99.504 


99.505 


V = 0.5 




TABLE 3.3 

COMPARISON OF NATURAL FREQUENCIES FOR A SIMPLY-SUPPORTED 

ORTHOTROPIC'*’ SQUARE PLATE 


Mode 

Ref <25> 

Present Method 

1 

0.0497 

0.04966 

2 

0.1120 

0.11205 

? 

0.1394 

0.15544 

4 

0.1987 

0.19874 


LT 


+ Ej_/Ej = 1.9, Glt/^T = 


e.44 



TABLE 3.4 


COMPARISON OF NATURAL FREQUENCIES OF A SIMPLY-SUPPORTED 
LAMINATE OF MATERIAL "A" <45* /-45* /45* /-45* ) . FOR 
VARIOUS ASPECT RATIOS 


Aspect Ratio 

12 

F.E.M 

Present Method 

0.2 

9.816 

9.556 

0.6 

15.689 

• 

CD 

1.0 

24.545 

24.109 

1.2 

29.521 

29.102 

1.6 

40.655 

40.414 

2.0 

55.989 

55.688 


a/h = 50. 




TABLE 3.5 

COMPARISION OF NATURAL FREQUENCIES OF CLAMPED SQUARE TAPERED 
ISOTROPIC PLATES FOR VARIOUS TAPER PARAMETERS tr. 

Taper Parameter, tr Ref. 24 Present Method 


0.2 

9.15 

9.442 

0.6 

10.185 

10.507 


0.8 


10.927 


10.728 





hb. of terms (noft) 


Fig. 3.1 : Convergence pattern for a simply supported 
anti-symmetric angle ply laminate; 
45/-45/45/-45; o/h = 50; 
aspect ratio = 1 . 






ratio o/n 

Fiq.3-2 : Variation of non-dimensionoi frequency with 
side/thickness ratio for simpi)/ supported 
anti -symmetric angle-ply laminates of 
material B; 45/ -45/ 4-5/ -45 
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Vc nation of non-dimensional frequency with 
aspect ratio for symmetric cross-ply 
!ominates;0/9G/90/G;(S.S.S.S);Q/h=100. 




Fig.3.4 ; Variation of non-dimensional frequency with 
aspect ratio for symmetric cross-piy 
ia nn inates; 0/ 90/9 D/ 0 ;(C. C.C.C);a/ h - 1 00. 
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Fiq3 5 : Variation of non-dinn,ens!onQi frequency With 
aspect ratio for anti-symmetric cross-ply 

laminates;0/9G/0/9G:(S.S.S,S.);a/ h-100 



a^ct ratio, p 


Fiq.3.6 : VariGtion of non-dtmensionol frequency with 
aspect ratio for anti-symmetric cross- ply 

laminates;0/90/0/90;(G.C.CX);a/h=1 00. 



fl^ct ratio, p 


Fig.3.8 : Variotion of non-dimensional frequency with 
aspect ratio for anti-symmetric angle-ply 
lanninates;30/-30/30/-30;(S.S.S.S);a/h=50. 






Fig.3.9 


Variation of non-dimensional frequency with 
aspect ratio for onti-symmetric angle-ply 
laminotes;i5/-45/45/-45;(S.S.S.S);a/h-50. 
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Fig 3 10 ; Variation of non-dimensiuonal frequency with 
aspect ratio for anti-symmetric angle-ply 

laminates;! 5/- 1 5/15/- f 5;(C.C.C.C);a/h-50. 


a^ct ratio,p 


yaridion of non-dimensional frequency with 
aspect ratio for anti-symmetric angle-ply 
laminates;50/-50/30/-50;(C.C,C.C);a/h = 50. 




D mctenai A 
0 material B 


a^ct ratio ,p 


Variation of non-dimensional frequency with 
aspect ratio for anti- symmetric angle-ply 
laminates;4-5/-45/45/-4-5;(C.CX.C);a/h = 50. 




Fig.3.13 : Variation of non-dimensional frequency with 

theta for simply-supported sguare anti-symmetric 
ongle-piy laminates of material A. 




fiber orientation , e 

Fig. 3.1 4 : Variation of non-dimensional frequency with 

theta for simply-supported square anti-symmetnc 
angle-ply laminates of material B. 




fiber orientation ,0 

Fiq3 15 : Variation of non-dimensional frequency with 

theta for clomped square onti -symmetric ongle-ply 
laminates of material A. 
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16 ; Variation of non-dimensiona! frequency with 

theta for damped square anti-symmeiric angle-ply 
laminates of material B. 
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CONCLUSIONS AND SUGGESTIONS FOR FUTURE WORK 

4.11 Conclusions: 

The following can be conclusively said about the present 
technique : 

► The present formulation is amenable to easy and extremely 
rapid programming. 

» The choice of the starting terms for the series is very 

easily made and the automatic generation and simple Gaussian 

integration possible, make the technique an elegant one. 

• Convergence is rapid with a fewer number of terms, due to 
the orthogonal i sat ion and this technique is excellent for a 
Ritz type analysis. 

E4. 2] Suggestions for future work: 

• For a particular type of problem, the program can be 

optimised by studying the nature of boundary conditions and 
avoiding repetitions of the term generation and integration 

operations. 

• The number of Gauss points to be used can be regulated by 
the number of terms being used in the series. 

• The present technique can be further exploited for free 
vibration analysis of thick plates if shear deformation is 
included in the formulation. 
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APPENDIX 


Component matrices of the sub-matrices that make up 
iffness matrix: 


[ al] : [a^ . I (3° } 


[a6^ 3 

2 2,0 o , 

: a p 1 cx^ ot^ 3 

* [/5^ ] 

[al6j 

2 r r 1 o 1 

: a p [a^ ] - 

1(3° + [oi° 

T7 T?-* '• ? 

[ 

2 T 1 /!jO 1 

: a p [ ot^ J - 

[ /5° cx^ 3 

<■ “‘■uv’ 

2 I* o 1 

: a P [a^ J 

. [ a° ] 

[al6 1 
'■ uv^ 

* a £ cx^ (3^ j , 

L ( ^ T> •* 

[ a26 .] 

2 2, o rtO 

; a p [ ot ^ 

] . ] 

J * 1- ' 7J J 

[a2 ] 

2 2 , «o , 

: a p 3 

r 1 1 1 

'■ 7? 7? •* 

Ca6^] 

S a^ c 3 . 

f O O ’ll 

7J 7? ■' 

t a26^] 

S a^p c ^^3 


|bl3 : 

-a[«^^|3 *1/5° ^ 3 

Cb6j 

_ 2 , o .1 1 

i -2ap I 4>^ 1 


C 

Zr 1 1 

; -ap [ 3 



the 



[ hisj : } • [/5° ] + [ cx° ] - [ 7?^ 3 ] 

[ b26^3 : [«° ^ I 

i b23 : -a/^’ lft° <P° } - [o.^ 3 

[ b6^3 : - 2aju [/5^ ) . [a° 4>^ ] 

[ bl2^3 : -agu [/?° ] . [a^ 3 

[ bl6^3 : ~a [ft^ 4>^ } - [oi° 3 

[ b26j : j^2[ 4>^ 1 - [ 3 * Ift^ 4>^ } - [ o° 3 ] 

C dl3 : 14>^ <P^ J - 3 

[ d23 : 4>° 1 . [ 4>^ 3 

[ d63 : 4^^ 3 - 3 

[dJ23 3 + [^° 3 .i4>l 4>°^^ 

Idl63 S [^1 4>l 3 - 1 

i:d263 * ^ 3 - ] + I4>l 4>'^ 3 . l] 



